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Nuclear magnetic resonance (NMR) spectroscopy has approached the limit of single molecule 
sensitivity, however the spectral resolution is currently insufficient to obtain detailed information 
on chemical structure and molecular interactions. Here we demonstrate more than two orders of 
magnitude improvement in spectral resolution by performing correlation spectroscopy with shallow 
nitrogen-vacancy (NV) magnetic sensors in diamond. In principle, the resolution is sufficient to 
observe chemical shifts in ~1T magnetic fields, and is currently limited by molecular diffusion at 
the surface. We measure oil diffusion rates of D = 0.15 — 0.2nm^//is within (5nm)^ volumes at the 
diamond surface. 


The extension of nuclear magnetic resonance (NMR) 
spectroscopy to the nanoscale m has important im¬ 
plications for scientific research and technology. In par¬ 
ticular, the attainment of sensitivities capable of detect¬ 
ing individual nuclei HE] has opened the way for non¬ 
destructive analysis of single molecules. In this context 
nitrogen-vacancy (NV) centres in diamond are particu¬ 
larly interesting given they have demonstrated outstand¬ 
ing nanoscale magnetic field sensitivity through a series 
of foundational experiments [IHIS]. However, whilst the 
sensitivity of NV-based NMR allows for few spin detec¬ 
tion, the spectral resolution is below that required to 
observe chemical shifts, thereby severely restricting the 
determination of chemical structure. 

The NV center is a paramagnetic defect in diamond 
consisting of a nitrogen impurity and a neighbouring va¬ 
cancy. The spin triplet ground state possesses long co¬ 
herence time and can be detected optically via coupling 
to optical transitions. Measuring the energy splitting be¬ 
tween magnetic sublevels of a single NV center allows 
determination of the local magnetic field [7]. The sensi¬ 
tivity of the measurement is then defined by the precision 
that the energy levels can be determined, which depends 
upon the coherence time and associated linewidth of the 
spin transition M- A narrower linewidth not only al¬ 
lows smaller field strengths to be detected, via shifts in 
the line position, but also allows finely separated energy 
levels to be resolved. 

For NV-based sensing, closely spaced energy levels may 
arise from dipolar interaction between the NV electronic 
spin and nuclei in molecules a few nanometers away, or 
from the internal energy structure of such molecules. 
Spectral resolution below part per million levels (with 
respect to applied magnetic field) is vital to observing 
chemical shifts, quadrupole interactions and dipolar cou¬ 
plings which allow molecular structure to be derived. 
Likewise, the identification of hyperfine interactions in¬ 
duced by NV centers at few nanometers distance from nu¬ 


clear spins requires sub-kilohertz spectroscopy. Therefore 
prolongation of the NV sensing time beyond this limit is 
important to enable applications of nanoscale sensing to 
structural analysis, and is especially relevant for centers 
located close to the diamond surface. 

One approach to improve NV magnetic resonance spec¬ 
troscopy is to remove decoherence sources which act to 
limit the detection time, either using quantum control 
techniques (e.g. dynamical decoupling) or by physical 
treatment (e.g. sample cooling and purification). Re¬ 
cently nanoscale NMR with shallow NV centers was re¬ 
alized, with spectral linewidths limited by the NV de¬ 
phasing (T 2 ) time [im]. Dynamical decoupling proto¬ 
cols were implemented to extend the NV dephasing time, 
resulting in a commensurate increase in sensitivity and 
spectral resolution. Dynamic decoupling alters the phase 
evolution of the sensing qubit such that the effect of un¬ 
wanted noise is removed, whilst retaining sensitivity to 
signals at a particular frequency. Pulsed decoupling tech¬ 
niques (e.g. XY-N, CPMG-N) operate by filtering signals 
at harmonics of the pulsing period whereas continuous 
decoupling (e.g. spin-locking, Hartmann-Hahn polariza¬ 
tion transfer) tune the driving frequency of the qubit 
[Fig-B a)]. Limitations in driving strength and fidelity 
however, reduce the effectiveness of decoupling, and pro¬ 
vide imperfect filtering. For example spectral noise faster 
than the decoupling speed is not removed, leading to pro¬ 
posals based on quantum error correction pffil [T6] . For 
NV sensing this problem is particularly salient, as poor 
scaling of dephasing time with number of applied pulses 
has been observed for shallow NV centers HZllIHI. As a 
result spectral resolution on the order of 10 - 100 kHz 
remains the benchmark, which is more than an order of 
magnitude worse than the intrinsic NV relaxation (Ti) 
time. 

An alternative avenue towards high resolution spec¬ 
troscopy is to transfer the signal encoded in the NV 
phase to spin population and then rely upon the long Ti 
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FIG. 1. Nanoscale magnetic sensing based on coherent phase and spin population timescales, (a) Magnetometry protocol using 
coherent phase accumulation: the sensing qubit is driven/pulsed at a desired frequency with e.g. an XY8 sequence. Evolution 
around the Bloch sphere is altered by magnetic fields at the sensing frequency, leading to differences in phase accumulation and 
finally spin population upon readout. Spectral resolution depends on the T 2 time of the sensor. Protocol using spin population 
correlations: an initial phase accumulation is stored in the sensing qubits spin population. A second phase accumulation is 
correlated with the initial measurement by a second mapping to spin population. The final population is a coherent sum of two 
measurements, dependent on the phase evolution of the external field during time, t. (b) A single nitrogen vacancy spin qubit 
below the diamond surface is used to measure the NMR spectrum of protons in a nanoscale volume of oil molecules. A magnetic 
field is applied with a permanent magnet and microwaves are carried by a copper wire, (c) Due to the extremely near-field 
interaction between the sensing qubit and sample, the recorded NMR spectrum is highly dependent on the NV distance to the 
surface. Diffusional motion and interactions near the surface of the diamond affect the recorded spectrum. 


time of the NV center. The idea of extending resolution 
to timescales limited solely by spin lattice relaxation is 
closely related to stimulated echo techniques m and was 
recently demonstrated for coherently coupled spins in di¬ 
amond [201 mi- Although this approach does not enhance 
the sensitivity of the NV center which relies upon coher¬ 
ent phase accumulation, spectral resolution can, in prin¬ 
ciple, be extended to the Ti limit [Fig.j^a)]. During the 
first “preparation” sequence consisting of a multipulse 
(e.g. XY8) echo, the magnetisation of processing nuclear 
spins is mapped onto the NV electronic spin population. 
After an evolution time t, during which the phase of the 
nuclear magnetisation progresses, the subsequent nuclear 
phase is then correlated with the initial phase by a second 
“read” echo sequence which again maps to the NV spin 
population [22]. The resultant signal is a modulation at 


the frequency of nuclear precession. 

A schematic of the measurement is shown in Fig.[^b). 
We use shallow NV centers implanted 2-5 nm into a (100) 
diamond surface by 2.5 keV N+ implantation. The di¬ 
amond is 99.999% isotopically enriched, with ^10 
part per billion impurity content so that the NV spectral 
environment is dominated by magnetic species at the sur¬ 
face or outside of the diamond m- A 400 Gauss mag¬ 
netic field was applied along the NV axis ((111) diamond 
crystal axis), resulting in a proton Larmor frequency of 
1.7 MHz (^H gyromagnetic ratio, 7 h = 4.25 kHz/Gauss, 
note the use of real frequency, not angular frequency val¬ 
ues). The NMR spectrum of statistically polarized pro¬ 
tons in an oil layer placed on the diamond surface was 
then measured with an XY8 sequence [2|. A peak in the 
spectrum corresponding with the Larmor frequency 
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FIG. 2. Comparison between coherent phase spectroscopy 
and correlation spectroscopy, (a) Proton NMR spectra ob¬ 
tained with XY8-4 (magenta), and XY8-8 (green) dynami¬ 
cal decoupling, giving a spectral linewidth of 110 kHz and 
74 kHz respectively. The spectrum of an externally applied 
AC field, measured with correlation spectroscopy (blue) has a 
corresponding linewidth of 0.47 kHz (inset), (b) Plot of data 
obtained with correlation spectroscopy before Fourier trans¬ 
form. The 1.7 MHz AC field was under-sampled in order to 
observe the envelope of the correlation signal. The Ti decay of 
the NV center (Ti = 1.7 ms) is shown as the grey dotted-line. 


records the root-mean squared magnetic field amplitude 
of the protons, as seen by the NV center [FigSa)]. From 
the amplitude of the nuclear spin signal [22f7 we calcu¬ 
late the NV distance to the diamond surface as d = 5 nm, 
with a detection volume of (5nm)^ [Fig.j^c)]. 

The resolution of the spectra shown in Fig.|^a) is ini¬ 
tially determined by the NV dephasing time. For an 
XY8-4 sequence consisting of 32 decoupling pulses, the 
measurement filter function has a linewidth of 110 kHz 
due to the ^ 10 /is NV T 2 time. By increasing the num¬ 
ber of decoupling pulses the NV T 2 time could be in¬ 
creased, leading to a narrower NMR signal with 74 kHz 
linewidth [Fig.j^a)]. However the resolution of these ini¬ 
tial measurements remains poorer than the expected pro¬ 
ton linewidth (typically in the tens of kilohertz range for 
solid-state NMR and less for liquid-state NMR), meaning 
little information about chemical structure or sample re¬ 
laxation rates can be obtained. We therefore performed 
correlation spectroscopy out to millisecond timescales, al¬ 
lowing a significant improvement in spectral resolution. 

As a demonstration of the technique we measured an 


effective delta frequency signal from a high-stability sig¬ 
nal generator. An XY8 sequence followed by a variable 
time delay and a second XY8 sequence resulted in os¬ 
cillations of the NV population at the signal generator 
frequency [22]. In order to detect the correlation oscilla¬ 
tions over long evolution times, whilst limiting the num¬ 
ber of datapoints required, we applied an under-sampling 
protocol which recorded the envelope of the correlation 
signal. The sampling rate was chosen to be 17 kHz (i.e. a 
sample every 58.859 /is) between 2/j^/n and 2fL /(n — 1), 
with fn = 1.704 MHz, /l = 1.699 MHz and n = 201 to 
fulfil the Nyquist criterion [23|. The result is shown in 
Fig.j^b), from which we obtained the linewidth plotted 
in Fig.j^a) by Fourier transform. A resolution improve¬ 
ment of over two orders of magnitude to 470 ± 40 Hz can 
be seen, in comparison to coherent phase sensing mea¬ 
surements. 

To our knowledge, the spectral resolution demon¬ 
strated here is a record for magnetometry with shallow 
NV centers [2llll[I71[l8l|2il[25]. Not only does the resolu¬ 
tion outperform current nanoscale NMR techniques, but 
it is comparable to measurements performed on protected 
NV centers deep inside diamond EUHillTl. To examine 
whether the linewidth is indeed limited by spin lifetime, 
we independently measured the Ti time of the NV cen¬ 
ter. We find the Ti decay of 1.7 ms closely matches the 
damping of the correlation signal [Fig.j^b)]. The good 
agreement between Ti decay and the measured linewidth 
indicates that it is the timescale which information is 
stored in the NV spin population rather than drifts in 
the external magnetic field or experimental setup, that 
limit the resolution of these measurements. 

NMR spectroscopy with few hundred hertz resolution 
is capable of identifying molecular structure in moder¬ 
ate magnetic fields based on chemical shifts. In Fig.[^ 
we show the simulated ^H and NMR spectra for 
acetic acid and methyl formate with 470 Hz resolution, 
where we have neglected any intrinsic sample relaxation 
and diffusion. The two molecules have identical chem¬ 
ical makeup, differing only in structure, meaning they 
cannot be distinguished by techniques sensitive to pro¬ 
ton/carbon ratio. We show in 5 Tesla magnetic fields, 
the chemical shift arising from different functional groups 
can be resolved, which allows for unambiguous identifica¬ 
tion by comparison to NMR databases [28| . By compar¬ 
ison, coherent phase spectroscopy is unable to yield any 
structural resolution due to the poorer linewidth [Fig.[^ 
(magenta)]. The higher chemical shift experienced by 
carbon atoms mean that chemical shifts can be resolved 
in magnetic fields of 1 Tesla, making NMR a more 
attractive option for NV based spectroscopy [Fig. Kb, 
inset)]. 

Acetic acid and methyl formate were chosen as they 
contain carboxylic acid and ester functional groups which 
are present at the diamond surface. Such groups are tar¬ 
gets for diamond functionalisation and expected to play 
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FIG. 3. Chemical shift resolution with an NV sensor, (a) 
Simulated NMR spectra of acetic acid and methyl formate 
at 5T with 470 Hz resolution from correlation spectroscopy 
(blue) and 110 kHz resolution from phase spectroscopy (ma¬ 
genta). (b) Simulated NMR spectra at 5T with correla¬ 
tion and phase spectroscopy. (Inset) Identical NMR spectra 
as (b), but for IT magnetic field, resulting in factor of five 
poorer resolution. Note, the spectra only display the intrin¬ 
sic resolution of the NV sensor, sample relaxation/diffusion is 
not taken into account. 

a role in the charge state and decoherence of shallow NV 
centers. We note that the realisation of chemical shift 
spectroscopy is predicated on the ability to both readout 
the NV spin state and preserve long Ti times at high 
magnetic fields, both of which have been demonstrated 
recently [29l [30]. Additionally, the intrinsic sample re¬ 
laxation rates must be low enough to allow for high res¬ 
olution spectroscopy. 

In order to determine whether detection of chemical 
shifts is indeed readily achievable, and to investigate 
relaxation mechanisms at the diamond surface we per¬ 
formed correlation spectroscopy of protons at the dia¬ 
mond surface with another NV center of similar depth. 
Fig.|^ (a) shows the signal measured to 27 /is, with char¬ 
acteristic oscillations at the ^H Larmor frequency. Again, 
we applied an under-sampling protocol to record the en¬ 
velope of the correlation signal. The sampling rate was 
chosen to be 0.81 MHz between 2fii/n and 2fi^/{n — 1), 
with fii = 1.75 MHz, = 1.65 MHz and n = 5. Os¬ 
cillations continue beyond 100 /is, which although longer 
than could be investigated with XY8 decoupling, is much 
shorter than the T\ time of this NV center [Fig.[^ (b)]. 
From the Fourier transform of the signal we determine 



FIG. 4. (a) Nanoscale proton correlation spectroscopy, (a) 
Oscillations due to Larmor precession of protons near the NV 
sensor, shown at two time intervals, (b) Long-time correla¬ 
tion signal obtained by under-sampling. The Ti decay of the 
NV center (Ti = 2.1ms) is shown as the grey dotted-line. 
Simulations comparing the decay rate for an oil diffusivity 
of 0.1 (magenta), 0.2 (blue), and 0.5nm^//is (green) with an 
oil density of 1.09 g/mL and a proton nuclear spin density 
of 50nm“^. (c) Fourier transform of the correlation signal, 
showing the measured proton NMR linewidth of 20 kHz (blue) 
compared to the intrinsic linewidth of the NV sensor (grey 
dotted-line), and conventional NMR proton spectrum of oil 
(green), (d) Linewidth broadening vs depth for 6 NV centers 
(red datapoints), and fit to data with D = 0.15 nw?/jis (blue 
line), measured with XY8 spectroscopy. Minimum detectable 
broadening (5 kHz) shown as dotted line. 


the proton linewidth to be 20 kHz, which is now no longer 
limited by the NV dephasing time [Fig. |4](c)]. 

Interestingly the NV measured proton spectrum is sig¬ 
nificantly broader than the ^ 40 Hz width spectrum mea¬ 
sured by conventional NMR [Fig.j^ (c)] (see [22] for de¬ 
tails). In nanoscale NMR, molecules can leave the detec¬ 
tion volume by molecular diffusion which reduces inter¬ 
action time with the sensor and results in a broadened 
linewidth. In contrast, molecular diffusion in conven¬ 
tional NMR leads to narrower linewidths through mo¬ 
tional narrowing. To test whether diffusion is responsi¬ 
ble for the linewidth broadening, we measured the pro¬ 
ton linewidth with six NV centers at different depths 
using XY8 spectroscopy and recorded broadening (be¬ 
yond the measurement resolution) as a function of sen¬ 
sor depth [Fig.|^(d)]. Despite the lower resolution asso¬ 
ciated with using XY8 spectroscopy, we observe a clear 
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increase in spectral linewidth as the NV depth reduces, 
which is fitted by a pure diffusion model (Acj = 2Dld?^ 
where Acj is broadening) to give a diffusion constant, 
D = 0.15 ± 0.04nm^//is. We note that our model does 
not take into account nuclear-nuclear interactions, inter¬ 
actions with spins at the diamond surface, or a static 
boundary layer of non-diffusing liquid which would be 
expected to give a different depth scaling. 

To validate the diffusional model we numerically sim¬ 
ulated the magnetic field produced by diffusing protons 
at the diamond surface. We used a cluster correlation 
expansion to determine the magnetic field produced by 
3000 randomly moving nuclear spin pairs, with their mo¬ 
tion calculated in 0.1/is timesteps [22]. As shown in 
Fig. 0 (b), a diffusion coefficient of 0.19 ± 0.05nm^//is, 
gives a best fit to the experimental results. Calculation 
based on pure diffusion, x‘^ 2Dt, where x‘^ is the mean 
squared displacement, reproduces this value well for dif¬ 
fusion through a (5nm)^ sensing volume, with a decay 
timescale of 65 /is, in agreement with Fig. 0(b). 

In summary, we have demonstrated that correlation 
spectroscopy can be applied to nanoscale NMR for im¬ 
proving the spectral resolution of diamond magnetome- 
try beyond values limited by the coherence time of elec¬ 
tron spins. The protocol is limited by the longitudinal 
relaxation time of the NV spin (several milliseconds at 
room temperature) leading to sub-kilohertz resolution, 
and in principle, allowing chemical shifts to be observed 
in high magnetic field experiments. Importantly, by stor¬ 
ing information in spin population, susceptibility to low 
frequency noise is reduced, which has been shown to be 
dominant at the diamond surface HaUHlEs]. In contrast 
to related relaxation spectroscopy techniques [HlSI], the 
linewidth does not depend on the NV T 2 time. 

We have also shown how molecular motion at the dia¬ 
mond surface affects the recorded spectrum, which can be 
used to determine molecular diffusion rates. In addition 
we highlight the importance of chemical attachment in 
order to achieve high resolution spectra. Improvements 
in stabilisation of molecules to the diamond surface, cou¬ 
pled with correlation spectroscopy at high fields, promise 
to yield information on molecular structure, morphology 
and dynamics at the single molecule level. 
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